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Abstract 

A new model of a Quantum Automaton (QA), working with qubits is proposed. 
The quantum states of the automaton can be pure or mixed and are represented by 
density operators. This is the appropriated approach to deal with measurements 
q-i and dechorence.. 

The linearity of a QA and of the partial trace super-operator, combined with 
the properties of invariant subspaces under unitary transformations, are used to 
minimize the dimension of the automaton and, consequently, the number of its 
working qubits. 

The results here developed are valid wether the state set of the QA is finite or 
not. 
t — ■ There are two main results in this paper: 1) We show that the dimension re- 

duction is possible whenever the unitary transformations, associated to each letter 
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of the input alphabet, obey a set of conditions. 2) We develop an algorithm to find 
out the equivalent minimal QA and prove that its complexity is polynomial in its 
dimension and in the size of the input alphabet. 
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1 Introduction 

In order to understand computation in a quantum context, it might be useful to introduce 
as many concepts as possible from the classical computation theory to the quantum case. 

One of these basic concepts concerns the functioning of finite automata in classical 
systems. To extend these concepts from classical to quantum systems, several models 
have been proposed [1, 2, 4, 5], most of them dealing with language recognition. This 
kind of automata is usually referred to as acceptors, in the sense that they are designed 
to indicate, by giving a yes or no output, whether a given input sequence does or does 
not possess the proper characteristics. 

In this article we propose a new model for a quantum automaton which acts as a 
transducer [6]. This QA is prepared in a given initial quantum state, then a classical 
device reads an input string of letters and apply the corresponding transformations to the 
quantum part. After all the letters in the string are read, the QA is observed (measured). 
Finally, the output of the quantum measurement is written in an exterior classical "tape". 
In this sense, our QA is similar to the measured-once Quantum Finite Automata (QFA) 
model introduced by [1], where the measurements are made only after a sequence of letters 
is read, and it differs from the measured-many QFA model introduced by [2, 7] that is 
observed after reading each letter. Since the outputs of this QA are probabilities, it can 
be seen as some kind of extension of the probabilistic transducers [8, 9]. Let us stress 
that, our QA does not have necessarily a finite state set and it is not used for purposes of 
language recognition. One important application of transducers is in the measurement, 
control and identification of systems. 

The study of such problems starts with the minimal-state equivalent automaton which 
suggests us that an important issue to take into account, when a given QA is defined, is 
how to minimize it. As far as we know such a question was never addressed in the context 
of QA. 

The main purpose of this article is to study the minimization procedure of a QA which 
differs from the procedure used with classical automata. 



The physical support of the QA that we are proposing in this paper is a system of n 
two-state quantum particles which encode information in the form of quantum bits, qubits. 
Each qubit is a unit vector of a 2-dimensional Hilbert space. Beside the formal reasons 
mentioned above, the minimization of a QA is an important question to be addressed 
because qubits are a very expensive resource. 

The minimization problem can be studied in two complementary ways. The first 
one, concerns the minimization of the dimension of the underlying Hilbert space, which 
consists in looking for an equivalent automaton working with the smallest number of 
qubits. Once the dimension minimization is achieved, another question can be asked, 
namely the possibility of minimizing the cardinality of the QA state set. 

Most authors represent the states of a quantum automaton by pure states. In general, 
a quantum system is not in a pure state. This may be due to the fact that we have only 
partial information about the system, or that the system is not isolated from the rest of 
the Universe, so it is not in a well defined pure state. In such cases we say that the system 
is in a mixed state which is represented by a density operator. The density operator is the 
most general and powerful way of expressing the state of a quantum system, namely when 
quantum measurements are performed [10, 11, 12], therefore we adopt this formalism to 
describe the states of a quantum automata. 

Another advantage of describing the states of a QA by density operators, is to deal 
with quantum measurements performed in subsystems of a larger quantum system. The 
approach developed in this paper to minimize the dimension of the QA uses the operation 
of tracing out that transforms density operators of a Hilbert space in density operators 
of another Hilbert space of smaller dimension. The transformed density operators are 
named reduced density operators . 

A QA is a special case of automata, this is, a QA is a linear automaton since its 
quantum states are vectors of a Hilbert space and the transition and output maps are 
linear transformations. Based on this linearity, we apply some known results of the theory 
of invariant subspaces to derive the necessary and sufficient conditions that the unitary 
evolution operators must obey, in order to minimize its dimension. 



This article is organized as follows: in section 2, the physics background is presented 
and the physical notation, used in quantum systems, is introduced. In section 3, the 
definition of our quantum automaton is presented and its functioning is explained. The 
conditions for reduction of its dimension are derived in section 4. The minimization algo- 
rithm is developed and its complexity is computed in section 5, . Finally, the conclusions 
are presented in section 6. 

2 Physics background 

A quantum physical system in a pure state is described by a unit vector in a Hilbert space, 
i.e., a vector space with a inner product. In the Dirac notation, the pure state is denoted 
by | ip). The quantum automaton, that we are considering, works with n-qubits which 
can be physically realized by n two-state particles. The Hilbert space of such a system 
is H.2 = C^ ' 1 ^", this is a 2 n dimensional complex vector space. Ti,^ is a tensor product 
of n Hilbert spaces, each one associated to a qubit: Tl% = 0" ^(j) where 'H2U) is the 
2-dimensional Hilbert space of qubit j. For each 'H2H) we choose a special basis set, called 
computational basis, consisting of two orthonormal states | ij), (ij = 0, 1). A basis set for 
H.2 is {<8>j=i I ij)jij — 0, 1}. A general pure state of n-qubits is a vector superposition 
of the type: | *) = Ei J ,..., in=0 ^i....i™®?=i \ ij), where £-„...^=o | W...*„ ?= 1- The 
transposed-complex conjugate of | \I/) is denoted by (\P |. The inner product between 
I ^) and I $) is denoted by (<& | \J/) = Yll u ..., in=0 V'n,...,^^*!,...,^- We define | ^)($ | 
to be the linear operator from Ti^ — > Ti^, known as the outer product of two vectors 
of 7^2) whose matricial representation, in the above mentioned computational basis, is 

®? = i®2=i<<il*X*lh*>- 

According to the postulates of quantum mechanics, the operations that we can per- 
form in a quantum system are represented by completely positive maps [10, 14, 15]. 
To our QA we are going to apply: a) unitary transformations, represented by unitary 
operators U acting in 7i r 2 l : U :| ^) — >| ^)' = U | $). The unitary transformations 
preserve the norm of the vectors in Ti 1 ^ ■ b) measurements of a given physical quan- 



tity A [13]. Such a quantity is represented by an hermitian operator A in Ti^, named 
observable quantity. The measurements can be performed in one or more qubits. The 
possible outcomes of these measurements are given by the eigenvalues a k (always real) 
of the operator A. The probability of the outcome at is Vk = Tr{P k \ ty)($? |}, where 
Pfc = X)f=i | al)( a k I) * s the projection operator in the gu degenerated eigenspace Sk, as- 
sociated to the eigenvalue Ofe. The set of vectors, {| a k )}, constitute an orthonormal basis 
set in Ti.2, obeying the eigenvalue equation: A | a\) = ak \ a k ); < k < n;0 < i < g k , and 
Tr stands for the trace of an operator. After a measurement, the quantum state of the 
system changes according to the outcome at becoming, | \& ) = P 



fc- 



l*> 



We have already mentioned that the most general state of a quantum system is not 
a pure state but rather a mixed state. We say that the system is in a mixed state, and 
assign with the system a probability distribution {p r ,\ ?Pr),12rPr = 1}, meaning that 
the system is in the pure state | ip r ), with the probability p r . A straightforward way of 
describing such a state is by using density operators p defined by [13], p = J2 r Pr I Vv)(Vv I- 
The density operators are: a. Linear operators of Ti^ — > Ti^', b. Hermitian, p = p<; c. 
Semi definite positive, i.e., for any vector | \P) G Ti^, (^ I P I ^) — 0; d. The trace, 
Tr{p} = E fc Ei(4 | p | 4) = 1. 

The set Cn n , of all linear operators of Ti^ — ► Ti^ is a vector space. Given any two 
vectors A and B of C-y^ we can define a inner product function by (A, B) = TrjA^B}, 
called the Hilbert- Schmidt inner product. With this inner product function, the vector 
space £n™, becomes a Hilbert space. To any unitary operator, U :| \&) — >\ ^f)' = U | ^/) 
on Ti.2, we can assign another unitary operator, tj : A — > A = U(A) = UAU^ on Cyi n - 
Operators like U, that transform operators into other operators are called super- operators. 

An important application of the density operators formalism is as a descriptive tool for 
subsystems of a composite quantum system. In fact, most physical systems are constituted 
by two or more parts (subsystems). Let us assume that the observable quantity A(l) 
measures only part (1), of a system, for instance the part composed by the n\ first qubits 
and call part (2) to the remaining n<i qubits, this is Ti^ = Ti.^ 1 <8> Ti-2 2 - 

The observable quantities A(l) on a subsystem (1) can be extended to a composite 
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system (1) + (2) as follows [13]: A(l) = A(l)® 1(2), where A(l) denotes the correspond- 
ing observable for the same physical measurement, performed on the composite system 
and 1(2) is the identity operator in Ti^ 2 ■ The action of this operator on the vector, 
| #(l))<g> | tf(2)) is defined by, 

A(l) [| tt(l)>® | *(2)>] = [A(l) | tt(l)>] ® [1(2) | tt(2)>] (1) 

The spectrum of A(l) in Ti^ is the same as the spectrum of A(l) in Ti.% 1 , but with 
all the eigenvalues degenerated in ?{%, even if none of them is degenerated in Ti^ 1 - The 
probability Vk of obtaining the outcome a^, when the observable A(l) is measured over 
part (1) of the system is given by Vk = Tr{P k (l)p}, where P k (l) = Pfc(l) ® 1(2) and 
Pfe(l) is the projector operator on the eigenspace £&. 

A natural question we can ask is, would it be possible to define a density operator for 
subsystem (1), that gives the correct probabilities Vk for the outcomes of the observable 
A(l)? 

It can be shown [14, 15] that there is a unique transformation of Cy^ into £ n n i which 
gives rise to the correct description of observable quantities for subsystem (1). This 
transformation, called trace out or partial trace and denoted Tr2{.}, is a completely 
positive map of C H n w «i which is the space of all linear operators: £-#« — ► £ w «i . It is 
defined for any pair of finite dimensional Hilbert spaces C-y^ and £ w «i , with n\ < n. 
The image of the partial trace of a density operator p e C-n n is denoted by p 1 = Tr 2 {p} 
and is named reduced density operator for subsystem (1), whose matrix elements are 
P l {hj) = 12k=ni+i{k\p\k) = Efc=„ 1+ iP(^! jk). It means that we are averaging over H^ 2 - 
Any quantum operator that does not operate on Ti^ 2 commutes with the partial trace. 
The reduced density operators have the same properties as the density operators. Let 
(A) to be a subsystem of ha qubits, containing the first m qubits and let (B) to be the 
remaining tib qubits, then ha = n — ub- The partial trace has the following property: 

V k = Tr{P k (l)p} = Tr A {P k (l)p A } = Tr 1 {P k (l)p 1 }. (2) 

Where p A = Trsip} is the reduced density operator of subsystem (A) and Trs{-} is the 
partial trace over subsystem (B). This equality shows that it is enough to know the reduced 



density operator p 1 (or p A ) in order to compute the probabilities of all the outcomes of a 
measurement over subsystem (1). 

3 A quantum automaton 

We say that a quantum system of n qubits works as a quantum automaton M. = 
(Ti%, po, E, U, Q, A(l), Q) when the following conditions are fulfilled: 

1. H2 is the underlying Hilbert space of dimension 2 n . 

2. po G C-)-c n is the density operator of the initial quantum state of the n qubits. 

3. E = {a} is a finite set of input symbols (the input alphabet). 

4- Q is the set of reachable states given by: Q = {p w : w G S*} ; with p w = XJ w (po) = 
U^poUjj, and where the family {V w } w eT,* of unitary operators ofTi^ is built as follows: 
i) U = {Uo-jo-gS; is a finite set of unitary operators indexed upon a G E; ii) U e = l(n); 

and Hi) \J wa = U ff U,„. 

5. d:ExQ — > Q is the transition map defined by 

$(&, Pw) = Pwa = Uff(P«,) (3) 

6. A(l) is the observable quantity to be measured. The possible outcomes of a 
measurement of A(l) are its eigenvalues a^: SpecA(l) = {a& : A(l)(|a^) ® 1(2)) = 
a k{Wk) ® 1(2))}- The diagonal representation of the observable is A(l) = Z)fc°fcPk(l)> 
where Pk(l) — J2f k =i I a k k )( a% k I ®1(2) is the projector into the eigenspace Sk of A(l). 

7. Q is a set of output symbols (the output alphabet) whose elements are defined by 
the following output map: 

8. u:Q^Q 

ou(w,a k )=Tr{P k (l)p w } (4) 

For each w G E* ; the set Q w = {uj(w,a k ) : a k G SpecA(l)} is the unique probability 
measure for the state p w and the output alphabet, Q = {Q w ,w G E*; J^k^i 10 ? a k) — 1}; 
is a set of probability measures. 



The physical performance of this QA is as follows: a classical exterior device can read 
the symbols a G E stored in an external classical "tape" . Then it applies the corresponding 
gate, represented by the unitary operator U CT , to the quantum state p w G C-n n -, °f the n 
two-state particles where the amount of information is stored. The dimension of the 
QA A4, denoted by dim[Ai], is the dimension of its quantum states p w : dim[Ai] = 
2 n x 2™ = N. After each input string of lenght w, the observable A(l) is measured and 
the corresponding output is given by the probability u(w,a,k) = Tr{Pk(l)p«,}, which is 
stored in an exterior classical device. 

The concepts of equivalent states and of equivalent automata can also be defined for a 
QA as follows. 

Two states p Wi and p w . of the state set Q of a QA At, are said to be equivalent, and 
we write p Wi ~ p Wj iff Lu(XJ w (p Wt ), a k ) = uj(tJ w (p Wj ), a k ) for every w G £*, More precisely, 
two states p Wi ,p W] G Q are equivalent, iff for every a k G SpecA.il) and for every w G E*: 

Tr{p Wi P k (l)} = Tr{p Wj P k (l)} (5) 

Tr{U w ( Pu ,jP k (l)} = Tr{U w (p„, 3 )P k (l)}. (6) 

We define a quantum behavior of a QA A4 over the signature (E, Q), as a map (3m{w) '■ 
E* — ► Jl(p w ) for every w G E*. Furthermore, two QA Ai and A4.4, over the same signature 
are said to be equivalent and we write M. ~ At a, when they have the same observable 
quantity A(l) and the same behavior. The equivalence between quantum automata is 
a bit more restrictive than between classical ones. The reason is because the output 
alphabet of a QA is determined by the probabilities of the outcomes of the measurements 
of an observable quantity A(l), therefore all equivalent QA must have in common this 
observable. In the present case, where the QA works with qubits, we cannot replace it 
by another QA with the same signature, if it doesn't work with qubits. Moreover, there 
is no QA equivalent to A-'f with less than n\ qubits, since the observable quantity must 
always be A(l). 

The fact that the dimension of the underlying Hilbert space of a QA is finite doesn't 
imply that the cardinality, or size, of its state set Q is finite. When the cardinal of Q is 



finite we call the QA a Quantum Finite Automaton. In such a case the output alphabet 
is also finite. 

The unitary operators U CT should have some properties in order that a QA might be 
also a QFA. 

Proposition 1: The number of states p\ = U'poU^ (I G N), generated by a unitary 
operator U ofH^, is finite iff there is a p & N and p < 2 n x 2 n , such that XJ P — 1. In 
this case the unitary operator U has p different imaginary eigenvalues Uj = e i27V ^ p , with 
1 < j < P, and the set generated by U is a finite group. 

Proof : The finite number of states of the QA is the result of the well known properties 
of any group of cyclic unitary operators. Q.E.D. 

Corollary 1: Let M. be a QA. If there is at least a a G E such that the set generated 
by the unitary operator U CT is not a finite group then the QA M. is not finite. 

Corollary 2: If for each a G £ the set generated by U CT is finite and if: \J ff i\J ff = 
Uo-Uo-/ for every a, o' G X, then {U w } we s* is a finite commutative group, and the QA Ai 
is finite. 

4 Dimension reduction 

Let us consider a QA, M. = (Ti.2,Po, S, U, Q, A(l), Q), composed by two subsystems (1) 
and (2) working with n\ and n2 qubits, respectively. The total number of qubits of the QA 
Ai is n = n\ + n 2 . The observable quantity of the QA is A(l) = A(l) ® 1(2) . Measuring 
this physical quantity corresponds to "inquiring" only about part (1) of the total system. 
A natural question we can ask is: In which conditions would it be possible to replace the 
QA M. by another equivalent QA M. a, using only ua {n\ < ua < n ), qubits? Such a 
QA would obviously have a smaller dimension than the original one. Moreover, if ua = n 
is the least number of qubits such that M.a — M- is equivalent to M. then, M. is the 
automaton with the minimal dimension or saying in other words, M. is the QA working 
with the minimal number of qubits. 
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What is usually understood by the minimization process of a classical automaton is 
the minimization of the cardinality of its sate set. However, this question can only be 
addressed if the QA is also a QFA. Instead, we are going to show that given a QA Ai, 
with a finite or infinite number of states, it is always possible to determine whether its 
dimension is minimal. 

A QA is a special case of automata, this is, a QA is a linear automaton since its states 
are vectors in a Hilbert space and the transition and output maps are linear transforma- 
tions. One of the distinct advantages of linear automata is that their behavior can be 
studied by analytical rather than the ennumerative techniques used in general automata. 

The procedure to minimize the cardinality of a QFA follows the usual approach de- 
veloped for classical automata, which consists in looking for equivalent states and the 
corresponding equivalence classes, however there are some specificities for quantum sys- 
tems that we will present in a future work. 

4.1 Disentangled subsystems 

It was mentioned in section 2, that the unique operation on a system that gives rise to 
the correct description of an observable quantity referred to one of its subsystems, is the 
partial trace over the remaining subsystems. This property will be the keystone in the 
process of dimension reduction of a QA. 

Theorem 1: Let M. = (Ti.%,Po, S,ZY, Q, A(l), $1) be a QA whose initial state po can 
be expressed in the form, p = p$ ® p$ , where p$ is the reduced density operator of part 
(A) of the system, working with ha qubits and containing subsystem (1), and p$ is the re- 
duced density operator of the remaining part (B) of the system, working with ub = n — n^ 
qubits. If all unitary transformations XJ a G U, had the form XJ a = XJ^~ ® U^ (where 
U^- and U^ are unitary transformations in £ n n A and C h ^b , respectively), then the QA 
•M-a — CH-2 A iPo j ^Ma, Qa, A(l), Q) is equivalent to M., and 

1 - H,2 A is the underlying Hilbert space with the minimal dimension 2 nA . 

2 - Pq = Trs{po} is the initial state. 

3 - The input alphabet, E = {o~}, is the same as in QA M.. 
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4- Qa = {pi-w G £*} with pi = Tr B {p w } = U^(p^) = U^/tfU^, and {U^}„, eS . is a 

family of unitary operators ofWf built as follows: i) XJf = 1(A); ii) Ua = {U^} CTe £ ; is 

a finite set of unitary operators indexed upon a G £ and Hi) U^, = U^U^. 

5 - 5 A : E x Q A -»• g A snc/i </io< p^ = ^(a, p^) = U^(p^) 

£ - The output function u>a '■ Q — ► ^ zs defined by UA{pi,a k ) = XVa{Pa;(1)p^} = 

uj(p w ,a k ). 

7-uj a : Qa—>& such thatuj A (pi,a k ) = Tr j4 {P k (l)p^} andfi = {fi w ;w G £*; Efc^A(w,a fc ) 

1} 7 to& Q w = {uja{Pw> a k) '■ &k G SpecA}, zs £/ie same as in QA M.. 

Proof: Using the properties of the partial trace operator: 
P A „ = Tr B {p a } = Tr B {tj a (p )} = Tr B {U^)U*(p B )} = tj^(p^). 
Using eq.(2), we obtain: u A (pi,a k ) = Tr A {Pfe(l)p^} = Tr{P fe (l)p„,} = u(p w ,a k ). 
Q.E.D. 

Theorem 1 says that if part (A) is initially disentangled from part (B) of the system 
and, if the unitary transformations preserve this disentanglement, then each subsystem 
(A) and (B) follows its own dynamics behaving as two separated systems. In this situation 
we can forget subsystem (B) and consider the QA restricted to subsystem (A), because 
M. ~ M.a- The equivalent QA M.^ is named a sober QA. Theorem 1 is a sufficient 
condition for dimension reduction. 

In the next sections we are going to use some known properties of invariant subspaces 
under unitary transformations, in order to derive the classes of equivalent quantum states. 

4.2 Linear transformations and invariant subspaces 

Let T : £-#>•> — f C n ^A to be the linear transformation defined by T = Tr B {.}. Let 
/C = ker(T) and X = im(T), to be the kernel and the image spaces of the transformation 
T. In Appendix A it is shown that, Im(T) = Cj^-a , therefore T is a linear transformation 
onto Cn,"A. Let Q C Cn n be the subspace orthogonal to /C, i.e. £^» = <20/C, then 

m-2 2 2 

dim[jCftri] = dim[Q] + dim[K], where dim[.} means dimension of a vector space. It is 
also true that dim[jC?{n] = dim[X] + dim\jC\. and we conclude that dim[Q] = dim[X}. 
Since the spaces Q and X have the same dimension they are isomorphic (I s4 theorem of 
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isomorphism) . 

Let us denote dim[C n ^\ = 2 n x 2 n = d, dim[Q] = 2 nA x 2 riA = q and dim[K] = 
2 n x2 n -2^ x2 nA = A; and let B Q = {Bf-J = l,...,g}, and B K = {Bf;j = q + 1, ..., d}, 
be two orthonormal basis sets for subspaces Q and /C, respectively. Their basis vectors 
obey the following orthonormal relations, 

(Bft B?) = (BfS B?) = * y ; (B?,B J K ) = (B«BJ») = (7) 

where (., .) is the Hilbert-Schmidt inner product defined in C-H n - 

The properties of any unitary transformation U are completely determined by its effect 
on the basis- B Q U B K '. Thus, if 

Z = X>B? + E ^Bf. (8) 

j=l j=l+q 

is any vector in £^n, the image of Z under U is 

C(Z) = J>tJ(Bj>) + E ^tJ(Bf) (9) 

therefore, if the images of the basis vectors are known, the image of any vector can be 
determined. The image of the basis- B Q U B , under U is, 

u(b?) = E^ Qb ?+ E 4 Kb ? ;j = i,-,? (io) 

i=l i=l+g 

u(Bf) = E4 QB ?+ £ 4* B * ;j = ? + i,-,d (ii) 

j=l j=l+(jr 



where, 



uf = (B?, U(Bj Q )) = Tr{B?fj(B?)} ;i = 1, ..., g; j = 1, ..., q (12) 

fig* = (B«, U(Bf)) = Tr{B?\J(Bf )} ; i = q + 1, ..., d ; j = q + 1, ..., d (13) 

«g X = (B^,U(B J ?))=Tr{Br t U(B J 9)} ;i = g + l,..,d ;j = l,...,? (14) 

fZ^ = (Bp,U(Bf))=Tr{Bp t U(Bf)} ;i = l, ...,</ ; j = ? + l,-,d (15) 

are the matrix elements of the super-operator U in the basis- B Q U B . 



13 
Replacing U(B- ) and U(Bj < ) given by (10) and (11) into eq.(9), we obtain 

u(z) = E(X>gS + E *$%}*? + E {E€% + E «r^> B f ( 16 ) 

i=i i=i i=i+9 j=i+9 i=i .7=1+3 

A possible relation between the subspaces Q and /C of £^n and the linear transformation 
U, is invariance. We say that Q (or /C) is invariant under U if, for every X G Q (or, for 
every Y G /C) implies U(X) G 72. (or U(Y) G /C ). We also say that <2 (or /C) is invariant 
under a set of linear transformations U„ (<r G £), if it is invariant under each member 
of the set. When two spaces, say Q and /C, such that Q0/C = Hu n i are both invariant 
under U, then we say that U is reduced (decomposed) by the pair (Q, /C). The following 
proposition [16] defines the structure of the matrix representing a unitary transformation 
U under which Q (or /C) is invariant. 

Proposition 2: Let C-u? = <2 © /C to be the direct sum of the subspaces Q and /C. If 
/C (or Q) is invariant under a unitary operator U , so is the complementary subspace Q 
(orJC). Then, the unitary transformation U is reduced by the pair (Q, /C) and the matrix 
representation of U in the basis-B Q U B is 

°=( U o Q tr*) (17) 

where the block diagonal matrices 

U QQ = [uf\i=,,..., q ; j= l,.., q (18) 

U = [Uij ]i=q+l,...,d;j=q+l,...,d (19) 

are unitary matrices, and the off' diagonal blocks are 

U QK = [u%% =q+l _ d ;J=1 ,..., 9 = (20) 

U KQ = [^ j Q } i =l,...,q; j =q + l,...,d = (21) 

Whenever U is reduced by the pair (<2,/C) then, U = U^^©!} 1 ^ 1 ^ is the direct sum 
of the unitary transformations U^^ and U KK defined on the subspaces Q and /C, res- 
pectively. The transformation U^^ describes what U does on Q and the transformation 
U KK describes what U does on /C. 
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Proposition 3: Let /C G Ch? t° be invariant under the unitary transformations XJ a 
and XJij. Then JC is also invariant under the transformation U a U&- 

Proof: The product of diagonal matrices is also a diagonal matrix. QED. 

What will the action of a unitary super-operator U on a general vector Z G Lu n be? 

Applying the linear super-operator T = Trs{-} to both sides of this eq.(8), and 
denoting by B] = T(Bj ); j = 1, ...q, we obtain 

The set B 1 = {B: 1 ; j = 1, ...,q} is an orthonormal basis set of £ n n A (Appendix A), and 
Z A = Tr B {Z} G £ h «a is the image of Z. 

The image of Z under U is Z = U(Z) given by eq.(9). Applying T to both sides of 
it, we end up with, 

Z A = v/2^ ± fe^gS + EL. ^i) B" ( 2 3) 



=9+1 

4 = 1 



where Z A = T(Z). 

We would like to know if there is a unitary super-operator \J A : C- n n A — > £ w ™a , such 
that for every state Z G £«»: 

tJ A T (Z) = TtJ (Z) (24) 

The action of any super-operator U on any Z A G C- n n A is, 

U A (Z A ) = v^^E^ [EL^il Bi 1 (25) 



where uf, = (B; 1 , U^B, 1 )) are the elements of tJ A relative to the £> 7 -basis. Replacing 
(25) and (23) in eq.(24) and recalling that B 1 is a set of linearly independent vectors, we 
obtain 

This is a system of q linear equations, where the unknowns are the q 2 elements u A { of 
matrix XJ A . There are (q 2 — q) linearly independent non-trivial solutions. We can choose 
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freely (g 2 — q) numbers u^ and the remaining q values being uniquely determined. Each 
of these solutions will depend on the components Zj of the vector Z, except when the 
coefficients of every Zj, in eq.(26), are simultaneously null, i.e., 

u A = uf i = l,...,q; j = l,...,q (27) 

u« Q = 0; t = l,...,q; j = q + l,...,d (28) 

Comparing these conditions with the statement of Proposition 2, we conclude that 
the invariance of subspace Q, under the action of the super-operator U, is a necessary 
condition for the existence of an operator XJ A G £ h »a, obeying eq.(24). In such case 
XJ A = U QQ , and the matrix tj QQ = [w§ <3 ]i=i,...,gy=i 1 ...,g is unitary. 

The following proposition summarizes these results. 

Proposition 4. Let /C C £%™ to be the kernel of the linear partial trace transformation 
T : C^™ — ► £ n ^A . /C is invariant under the unitary operator U : £^« — ► £%« iff 
-Q K Q — o ( or U ^ = 0) . Then U = U^ © U KK and there is a unitary transformation 
fjA = tjQQ : C h «a -»• C H n A , such that, for every ZeC n ^, U A f (Z) = TU(Z). 

4.3 Condition for dimension reduction 

Applying the transformation T = Trs{-} to each sate of the state set Q = {p w = U w (po) : 
to e E*} of the QA A4, we obtain the set Qa = {pt — Trsipw} '■ w G £*}, where p A 
are the reduced density operators in £ n n A. In particular, the initial state po is mapped 
into p A = Trs{po}- The set Qa will be a set of reachable states of a QA M.a, iff there is 
a family of operators Ua = {U A : w G E*}, such that p A = XJ^p A XJ^ . We have shown 
in Proposition 4 that when U^ Q = 0, then U CT = U^ Q ©U^ CK and there is a unitary 
operator U A = U? Q . By Proposition 3, U w = U^ Q ©U* K (w G E*) iff U^ Q = 0, for 
every a G E. 

Moreover, the quantum automata .M and .Ma will be equivalent, if they have the 
same behavior. The probability of obtaining an outcome a k when the QA Ai is in the 
state p w is: u(p w ,a k ) = Tr{P k (l)p w } = Tr A {Pk(l)pi} = uj A (p^,a k ), where we used 
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eq.(2). This last equality shows that the quantum automata Ai and A4a have the same 
behavior: Pm( w ) — Pm a ( w ) '■ ^* ~^ ^> f° r ever Y w G S*. 

The results that we have been deriving can be summarized in the following theorem. 

Theorem 2: A tuple Ma = {Ti-2 A ,Po, E, Ua, Qa,A(1),Q), w#i n^ < n and 
whose state set Qa = {pi = 7Vb{U to (po)} '■ w G E*} is a QA equivalent to the QA Ai, 
iffXJ^™ = for every unitary operator Uo-es- In such a case, XJ a = XjQQ $tJ^ K : 

1 - H,2 A is the underlying Hilbert space with dimension 2 nA < 2 n . 

2 - Pq = Ttb{po} G £-H nA ^ s th e initial state. 

3 - The input alphabet, E = {cr}, is the same as in QA M.. 

4- The state set Q A = {pi : w G E*} with pi = Tr B {p w } = \J^(p w ) = U A (p£) = 
XJipiXJi^ , where {U^}„, e s* is a family of unitary operators of H.2 A built as follows: i) 
Ua = {U^"} ctGS; is a finite set of unitary operators indexed upon a G E ; it) XJf = 1(A); 
andm)\Ji w = \Ji\Ji. 

5 - The transition map ^:Ex Q A — *■ Qa su ch that pi a = 5a(°~, pi) = U A (p^). 

6 - The output function uja '■ Q — > ^ is defined by u>A(pi,a,k) — T r A{Pk(^)pi} — 
u(p w ,a k ). 

7 - The output set, Q = {Q w ;w G E*; J2k u A(w, a k ) = 1}, with Q w = {u>a(w, a k ) : a k G 
SpecA}, is the same as in QA M.. 

This theorem tells us how to built the physical support of the equivalent automaton 
M.a- It is a system of n^^i < ua < n) two-state quantum particles prepared in a 
quantum state pi = Tr B {p }, and submitted to quantum gates represented by the set 
of unitary operators U A = {U^} CTeS such that U A (p A ) = U^p^U^ 1 " and tj A = tj^ Q . 
The probabilities of the possible outcomes of a measurement of the observable A(l) = 
A(l) <S> 1(B), performed on the first n\ of the ua two-state particles, give the output 
symbols of the smaller QA A4a, with dimension 2 nA x 2 nA . 

If there is at least one input a G E such that U^ ^ 0, then there is no QA with 
dimension 2 nA x 2 nA , equivalent to M.. 
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5 Minimization algorithm 

Theorem 2 says that a QA M.^, working with ua < n is equivalent to the initial QA 
Ai iff the kernel of the partial trace transformation is invariant under all the unitary 
operators U CT , with uGE. However, it does not say whether ha is the minimum number 
of qubits n. In this section we present an algorithm to determine the minimal QA. 

Usually, the density matrices p and the unitary operators U CT , are expressed in terms 
of the computational Dirac basis-B^(n) of £#«, (Appendix A), rather than in terms of the 
basis-£> s U B , that is needed to check the invariance of subspace /C under the operators 
U CT . Therefore, the first step of the algorithm consists in constructing two ortonormal 
basis sets for the subspaces Q and /C, as explained in Appendix A. 

We assume that the n\ first qubits of the QA Ai belong to subsystem (1) where the 
observable A(l) is going to be measured and that we have already checked if the conditions 
of Theorem 1 are verified. If so we obtain the corresponding sober automaton. 

Given a QA M. = (H^, po, S, U, Q, A(l), Q) over a finite signature (E, T), we compute 
the QA M. = (Ti,2,Po, S,W, Q, A(l), Q), working with the minimal number of qubits n, 
as follows: 

1. Replace each Dirac vector <S> r l l = iV i ->;ij,i j = 0, 1, of the Dirac computational basis 
set B^(n) by the 2 n x 2 n matrix B r (r = 2 n (l — 1) + /'), whose elements are all null 
except the element b t f , that is equal to 1 and obtain the ordered computational basis set 
B%(n) = {B r ,r = l,...,2 n x 2 n }. 

2. Do Ha = n\. 

3. Replace each Dirac vector <8)?=ifj.j' £ &d{ n A)'-,ij-,ij = 0, 1, of the Dirac compu- 
tational basis set B%(n A ), by the 2 nA x 2 nA the matrix Bj, (r = 2 nA (l - 1) + I'), whose 
elements are all null except the element b t f , that is equal to 1 and obtain the ordered 
computational basis set B^{riA) = {B^r = 1, ... 1 2\ x 2 1 \}. 

4. Compute T(B r ) for r = 1, ...,2 n x 2 n . If T(B r ) = add the vector B r to the set 
S . If T(B r ) = B] 1 , with j = 1, ...,2 nA x 2 nA , add B r to the set Sj. 



18 



5. With vectors of each set Sj = {B rk (j) : T(B ric (j)) = Bj 1 ; k = l,...,2 n ~ nA } , 

jQ = l 

j ' y/2 n - n A 



j = 1, ..., 2 UA x 2 nA , obtain the following linear combinations: Bj = , n 1 _ n Y^C=\ A B rk (j) 



and build the set B Q = {Bj 4 : j = 1, ..., 2™ A x 2™ A }. 

6. With the vectors of each set Sj, j = 1, ...,2 nA x 2™ A , obtain the following linear 
combinations: Bu = y)?_i A CfcB r . (?) where c& = — -, — = for k — 1, ...2 n_ ™ A — 1 and 

J — ^ - 1 - k w / \/ t 2 Tl ~' n A 1 

Cnn — nj^ — — -L . 

7. Compute the vectors By corresponding to all cyclic permutations of the coefficients 
c k (k > 1) in eq.(31), and build the 2 nA x 2 nA sets Cj = {By;/ = l,...,2 n ~ nA - 1} with 
j = l,...,2 nA x 2 nA . 

8. Apply the Gram-Schmidt algorithm to each set Cj in order to obtain a set of 
orthonormal vectors Cf . 

9. Apply an appropriate ordering algorithm to the vectors given by (32) to obtain an 
ordered set of orthonormal vectors B K = {B^; r = 2 nA x 2 nA + 1, ..., 2 n x 2 n } and build 
B c = B Q U B K . 

10. Write the transition matrix C from the basis-i3^*(n) to the basis-i^ . 

11. For each XJ' a E U compute u^a) = Tr{Bf U^BJ^U^} with i,j = 1, ..., 2 n x 2 n . 

12. Build the matrices \j' a = [uij(cr)], i,j = 1, ...,2 n x2 n ,aG £. 

13. For each a G £ compute U CT = C _1 tJ CT C 

14. For each a e S compute wg x (a) = Tr {B^U^B^U^}, i = 2™ A x2 nA + l, ..., 2 n x2 n 
and j = l,...,2 nA x 2™ A . 

15. If there is at least one Uij(ap K ^ 0, go to 17. If for every a G £ and for every 
i = l,...,2 nA x 2 nA ,j = 2 nA x 2 nA + l,...,2 n x 2 n , u% K {a) = 0, goto 16. 

16. For each a G E build the matrices U^ Q = [u^{p)\ = Tr{W (Bf)B^};i,j = 

l,...,2 nA x 2 nA , go to 20. 

17. Do riA = ua + 1- 

18. If n^ = n, go to 21. If not, go to 3. 
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19. The initial QA M., is the minimal one. 

20. Write n = n A . 

21. Compute pb = Tr^jpo}- 

22. Write the set U = {U^ Q ; a e E}. 

23. The QA .M = (H%,p , E,ZY, Q, A(l), fi), with n < n qubits, is the minimal 



one. 



The complexity of this algorithm is computed in terms of the dimension of the ini- 
tial QA Ai, and of the dimension of the QA M-a which are dim[M] = 2™ x 2 n = N and 
dim[M.A\ — 2™ A x 2™ A = N A respectively. Another variable playing a role in the computa- 
tion of the complexity of the algorithm, is the size |E| of the input alphabet, since it gives 
the number of unitary operators U CT (gates) to be computed in the B c = B Q U S -basis. 
The cardinality of the state set Q, doesn't play any role in the minimization algorithm. 
The dimensions of the observable A(l) and the size \Q\ of the output set, are irrelevant 
variables in the algorithm. 

Theorem 3: Given a QA M. over a signature (E, Q), the minimization algorithm 
requires 0(p(|E|, dim[A4})) arithmetic operations. 

Proof: Here is the detailed analysis of the algorithm (bounding for each step the 
worst case execution time in terms of arithmetic operations): 

1 - 0(N 2 ) since the cost of ordering iV objects is quadratic and the cost of computing 
the matrix representation of iV vectors ®" =1 v. ■> ; i~, %, = 0, 1 is iV . 

2 - 0(1) since this is the cost of a constant. 

3 - 0(N\) since the cost of ordering Na objects is quadratic and the cost of computing 
the matrix representation of N A vectors ®"=i v-j ; ij, i, = 0, 1 is N A ■ 

J j j J 

4 - 0(N 3 / 2 N A ~ ) since the body of the cycle is running O(N) times and each run 
computes the partial trace T{B r } with the cost 0((N/N A ) 1 / 2 ) 

5 - O(N) since this is the cost of a sum of vectors of dimension N. 

6 - O(N) for the same reason as in 5. 
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7 - 0((NN A ) 1/2 ) since 0((N/N A ) 1/2 ) is the cost of a permutation of (N/N A ) 1/2 objects 
and there are Na permutations to compute. 

8 - 0{N^ 2 N A /2 ) since 0{N 2 (N/N A ) 1/2 is the cost of applying the Gram-Schmidt 
process to a set of (N/Na) 1 ^ 2 vectors of dimension N each, and this process is repeated 
for each of the N A sets Cj. 

9 - 0((N — Na) 2 ) since the cost of the ordering algorithm is quadratic in the number 
of objects and this is the number of vectors in the set given by eq.(30). 

10 - 0(N 2 ) since each of the N vectors of the £>^-basis has N coefficients, being N 2 
the dimension of the super-operator C. 

11 - 0(N 5 ) since there are N 2 operations corresponding to two loops, being the cost 
of the body of each loop 0(N 3 ). The body of each loop computes the matrix product 
{B^U^B; Uj"} that has a cubic cost, followed by the computation of its trace that has 
cost O(N). 

12 - 0(|X|iV 5 ) since the body of the cycle is running 0(E|) times and the cost of each 
run is 0(N 5 ). 

13 - 0(|E|iV 6 ) since the body of the cycles is running 0(E|) times and the cost of each 
run is cubic (product of matrices) in the dimension of the super-operators C which is iV 2 . 

14 - 0(|E|(7V - N A )N 3 ) the same reason as in step 12. 

15 - 0(|E|iV 2 ) since the body of the cycles is running 0(|E|iV 2 ) times and the cost of 
each run is 0(1). 

16 - 0((\Tj\NaN 3 ) the same reason as in step 12. 

17 - 0(\T.\log(N/NA)N 6 )) since the body of the cycle is running 0(log(N/NA)) times 
and the cost of each run is 0(|S|iV 6 ). 

18, 19 and 20 - 0(1) since this is the cost of a constant. 

21 - 0(N 3 / 2 Na ), the same reason as in step 4. 

22 and 23 - 0(1) since this is the cost of a constant. 
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In conclusion: all steps of the minimization algorithm requires a polynomial number of 
arithmetic operations in |E|, N = dim[Ai], Na = dim[M.A\ and the dim[M.A\ < dim[Ai]. 
QED. 

6 Conclusion 

In this paper a new model for a QA working with qubits was proposed and the problem 
of minimizing its dimension was solved. 

The quantum states of the QA were represented by density operators which is a power- 
ful approach to deal with measurements concerning only part of a quantum system. The 
linearity of the automaton and of the partial trace super-operator were used to derive 
the conditions for dimension reduction. It was shown that this is possible, if the kernel 
of the partial trace transformation is invariant under each of the unitary transformations 
associated to each letter of the input alphabet. 

It was also developed a minimization algorithm and it was shown that its complexity 
is polynomial in the size of the input alphabet and in the dimension of the QA. 

Let us stress again that the minimization of the dimension o a QA and the min- 
imization of the cardinality of its state set are different issues. While the dimension 
minimization can be addressed whether the QA is finite or not, the cardinality minimiza- 
tion is only possible if the automaton is a QFA. This problem will be studied in a future 
work. 

Finally, we refer the possibility that the minimization technique here developed, can be 
adapted to other kind of automata, namely Quantum Cellular Automata that implement 
quantum computation on qubits using spins [17, 18]. 

Appendix A 

The computational basis set for C-^ is an orthonormal basis set. In the Dirac notation it 
is expressed by B%(ri) = {^>k=i v i k i' 'ih,h = 0' 1} w ^ n v i i' = l*fc)(*fcl an d dim[B%(n)] = 
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2 n x 2™. To order the Dirac vectors of this basis we associate to each Dirac vector 
®fc =1 v i j G Bd(n); ik, i^ = 0, 1, the 2 n x 2 n matrix B r , whose elements are all null except 
the element b u >, that is equal to 1 and where r = 2™(/ — 1) + /'. Applying the same 
reasoning to the Dirac vectors ^>k=i v i % e Bd( n A)',ik,ik = 0,1, (which is the Dirac 
computational basis for the image space C>-u nA ) we associate to each vector <S)?=ifj j G 

2 3 : j 

Bj(nA)',ik, ^ = 0, 1, a 2™ A x 2 nA matrix B*, whose elements are all null except the element 
6 M ', that is equal to 1 and where r = 2 nA (l - 1) + /'. Call B 1 = {B^r = 1, ...,2 nA x 2 nA } 
to the ordered computational basis-set of the image space £ n n A . 

The image of any vector of the basis- B^'(n) by the transformation T = Ttb{-} is, 

n \ ua n 

v SkS ' k ® v ak< =®v ; rr fl { <g> «.„;} (29) 

yfc=l £=71^4+1 y fe=l fc=n j 4 + l 

Computing the partial trace we obtain, 




if 3j : s fc ^ s' fe 



n 

V 
fc=l ™ k=n A +l "I k w « =i "s fc s fc 



Let us denote by S = {B rk : T(B rk ) = 0; k = 1, ..., 2™(2 n -2™ A )} the set of all vectors 
of B^'(n) transformed by T into zero. Let us denote by Sj = {B rk (j) : T(B rk (j)) = 
B: 1 ; k = 1, ..., 2™ _ ™ A } the set of the vectors of B^(n) with the same image B, 1 . There are 
2 n ~ nA vectors B rk (j) G B£(ra) and there are 2™ A x 2™ A sets S,-. 

A straightforward way of obtaining an orthonormal basis set for the subspace Q con- 
sists in taking the following linear combinations of all vectors B rk (j) G Sj: 

i Y a ~ n A 

B? = TJPS g B, t 0) (30) 

The set B Q = {B- : j = 1, ...,2™ A x 2 nA } is an ordered set of orthonormal vectors and 
constitute a basis set for the subspace Q. In fact, applying the super-operator T to both 
sides of eq.(30) the result is the vector Bj = \ / 2 n ~ nA ~B'*. 

With the vectors of the sets Sj and So it is possible to build up an orthonormal basis 
set for the kernel subspace /C, complementary to Q. It is simple algebra to prove that the 
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following linear combination of all vectors B rk (j) G S 



j 



2 n - n A 



B lj= £ C ^ B r k (j) (31) 

k=l 

with Cfc = 2 n-i 4 _ 1 ; fc = 1, ...2 n ~" A — 1 and c 2 ™-™^ = — 1, is applied on by the super- 
operator T. 

Making all the cyclic permutations of the last 2 n ~ nA — 1 coefficients C& of eq.(31) it is 
possible to generated 2 n ~ nA — 1 different vectors B^- with I = 1, ..., 2 n ~ nA — 1. The vectors 
of each set Cj = {B y ; I = 1, ..., T~ nA - 1}; (j = 1, ..., 2 n ° x 2 nA ), are orthogonal to the 
vectors of the basis- £> s and to the vectors of the set Sq. 

The vectors of each set Cj are not mutually orthogonal. They can be transformed in 
an set of orthormal vectors, Cj~, by applying to each set Cj the Gram-Schmidt algorithm. 

Giving all these properties of the vector sets So an d Cj~, it is now clear that the 
following union of sets 

S uf = Y 2nA Cf. (32) 

is an ortonormal basis set for the kernel subspace /C. Applying an appropriate ordering 
algorithm to the vectors given by (32) we obtain an ordered set of orthonormal vectors 
B K = {Bf; r = 2 nA x2 nA + l, ..., 2 n x 2™}. 

Let B c = {B r ; r = 1,...,2" x 2™} be the ordered orthonormal computational basis 
set in Cf{n. This basis set is related to the orthonormal basis set £> s U B K through the 
equations 

2 n x2 n 2 n x2™ 

BP = E ^B s ; Bf=Y, c£B s (33) 

S=l 8 = 1 

where c% (i = 1, ..., 2 nA x 2 nA ; s = 1, ..., 2™ x 2") and c^ (i = 2 nA x 2 nA + 1, ..., 2 n x 2 n ; s = 
1, ..., 2 n x 2 n ) are the elements of the transition matrix C from the basis- B c to the basis- 
B Q U B K '. The action of U relative to the i3 c -basis is 

2 n x2" 



U'(Bj) = E ««Bi (34) 



i=l 
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where u^ = (U (Bj),Bi) = TV{U t(Bj)Bi} are the elements of U relative to this basis 
set. If U is the matrix of the super-operator U relative to the B c -basis, then 

U = C^tJC (35) 

is the matrix of the same super-operator relative to the basis- £> s U B . The matrix U 
given by eq.(35), is the one that must be used to compute the matrices U^ and XJ KCi 
of eqs.(18) and (21). 
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